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Abstract. We consider the problem of constructing bounded-degree planar geometric 
spanners of Euclidean and unit-disk graphs. It is well known that the Delaunay subgraph 
is a planar geometric spanner with stretch factor del ~ 2.42; however, its degree may 
not be bounded. Our first result is a very simple linear time algorithm for constructing 
a subgraph of the Delaunay graph with stretch factor p = 1 + 2-7r(fccos f-) -1 and degree 
bounded by k, for any integer parameter k > 14. This result immediately implies an 
algorithm for constructing a planar geometric spanner of a Euclidean graph with stretch 
factor p ■ del and degree bounded by k, for any integer parameter k > 14. Moreover, the 
resulting spanner contains a Euclidean Minimum Spanning Tree (EMST) as a subgraph. 
Our second contribution lies in developing the structural results necessary to transfer our 
analysis and algorithm from Euclidean graphs to unit disk graphs, the usual model for 
wireless ad-hoc networks. We obtain a very simple distributed, strictly-localized algorithm 
that, given a unit disk graph embedded in the plane, constructs a geometric spanner with 
the above stretch factor and degree bound, and also containing an EMST as a subgraph. 
The obtained results dramatically improve the previous results in all aspects, as shown in 
the paper. 



Introduction 

Given a set of points P in the plane, the Euclidean graph E on P is defined to be 
the complete graph whose vertex-set is P. Each edge AB connecting points A and B is 
assumed to be embedded in the plane as the straight line segment AB; we define its cost 
to be the Euclidean distance |^4-E>|- We define the unit disk graph U to be the subgraph of 
E consisting of all edges AB with \AB\ < 1. 

Let G be a subgraph of E. The cost of a simple path A = Mn, Mi, M r = B in G is 
X^=o \MjMj+i\. Among all paths between A and B in G, a path with the smallest cost is 
defined to be a smallest cost path and we denote its cost as cg(A, B). A spanning subgraph 
H of G is said to be a geometric spanner of G if there is a constant p such that for every 
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two points A, B £ G we have: ch(A, B) < p- cg(A, B). The constant p is called the stretch 
factor of H (with respect to the underlying graph G). 

The problem of constructing geometric spanners of Euclidean graphs has recently re- 
ceived a lot of attention due to its applications in computational geometry, wireless com- 
puting, and computer graphics (see, for example, the recent book [13] for a survey on 
geometric spanners and their applications in networks). Dob kin et al. [9] showed that the 
Delaunay graph is a planar geometric spanner of the Euclidean graph with stretch factor 
(1 + v / 5)tt/2 ps 5.08. This ratio was improved by Keil et al [10] to Cdei = 2-7r/(3cos (vr/6)) ~ 
2.42, which currently stands as the best upper bound on the stretch factor of the Delaunay 
graph. Many researchers believe, however, that the lower bound of ir/2 shown in [7] is also 
an upper bound on the stretch factor of the Delaunay graph. While Delaunay graphs are 
good planar geometric spanners of Euclidean graphs, they may have unbounded degree. 

Other geometric (sparse) spanners were also proposed in the literature including the 
Yao graphs [16], the 0-graphs [10], and many others (see [13]). However, most of these 
proposed spanners either do not guarantee planarity, or do not guarantee bounded degree. 

Bose et al. [21 13] were the first to show how to extract a subgraph of the Delaunay 
graph that is a planar geometric spanner of the Euclidean graph with stretch factor ~ 10.02 
and degree bounded by 27. In the context of unit disk graphs, Li et al. [HI [12] gave a 
distributed algorithm that constructs a planar geometric spanner of a unit disk graph with 
stretch factor Cdei', however, the spanner constructed can have unbounded degree. Wang 
and Li |14[ [To] then showed how to construct a bounded-degree planar spanner of a unit disk 
graph with stretch factor rreax{7r/2, 1 + 7r sin (a/2)} • Cdei and degree bounded by 19 + 27r/a, 
where < a < 2-7r/3isa parameter. Very recently, Bose et. al [5] improved the earlier 
result in [21 [3] and showed how to construct a subgraph of the Delaunay graph that is a 
geometric spanner of the Euclidean graph with stretch factor: max{ir/2, 1+7T sin (a/2)} -Cdei 
if a < it/2 and (l + 2\/3+37r/2-(-7r sin (it /12))-Cdei when tt/2 < a < 2ir/3, and whose degree 
is bounded by 14 + 2n/a. Bose et al. then applied their construction to obtain a planar 
geometric spanner of a unit disk graph with stretch factor max{Tr/2, 1 + 7rsin (a/2)} ■ Cdei 
and degree bounded by 14 + 2-7r/a, for any < a < tt/3. This was the best bound on the 
stretch factor and the degree. 

We have two new results in this paper. We develop structural results about Delaunay 
graphs that allow us to present a very simple linear-time algorithm that, given a Delaunay 
graph, constructs a subgraph of the Delaunay graph with stretch factor l+2ir(k cos (ir/k))~ l 
(with respect to the Delaunay graph) and degree at most k, for any integer parameter 
k > 14. This result immediately implies an O(relgn) algorithm for constructing a planar 
geometric spanner of a Euclidean graph with stretch factor of (1 + 2n(k cos (ir/k))^ 1 ) ■ Cdei 
and degree at most k, for any integer parameter k > 14 (re is the number of vertices in the 
graph). We then translate our work to unit disk graphs and present our second result: a very 
simple and strictly-localized distributed algorithm that, given a unit-disk graph embedded in 
the plane, constructs a planar geometric spanner of the unit disk graph with stretch factor 
(1 + 27r(A:cos (Tr/k))^ 1 ) ■ Cdei an d degree bounded by k, for any integer parameter k > 14. 
This efficient distributed algorithm exchanges no more than 0(n) messages in total, and 
runs in O(AlgA) local time at a node of degree A. We show that both spanners include a 
Euclidean Minimum Spanning Tree as a subgraph. 

Both algorithms significantly improve previous results (described above) in terms of the 
stretch factor and the degree bound. To show this, we compare our results with previous 
results in more detail. For a degree bound k = 14, our result on Euclidean graphs imply 
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a bound of at most 3.54 on the stretch factor. As the degree bound k approaches oo, our 
bound on the stretch factor approaches C^ei ~ 2.42. The very recent results of Bose et 
al. [5] achieve a lowest degree bound of 17, and that corresponds to a bound on the stretch 
factor of at least 23. If Bose et al. [5] allow the degree bound to be arbitrarily large (i.e., to 
approach oo), their bound on the stretch factor approaches (n/2) ■ C^ei > 3.75. Our stretch 
factor and degree bounds for unit disk graphs are the same as our results for Euclidean 
graphs. The smallest degree bound derived by Bose et al. [5] is 20, and that corresponds to 
a stretch factor of at least 6.19. If Bose et al. [5] allow the degree bound to be arbitrarily 
large, then their bound on the stretch factor approaches (vr/2) • C^i > 3.75. On the other 
hand, the smallest degree bound derived in Wang et al. [HI [15] is 25, and that corresponds 
to a bound of 6.19 on the stretch factor. If Wang et al. [HI [15] allow the degree bound to 
be arbitrarily large, then their bound on the stretch factor approaches (n/2) ■ C^ei > 3.75. 
Therefore, even the worst bound of at most 3.54 on the stretch factor corresponding to our 
lowest bound on the degree k = 14, beats the best bound on the stretch factor of at least 
3.75 corresponding to arbitrarily large degree in both Bose et al. [5] and Wang et al. [T4"l [15] ! 

1. Definitions and Background 

We start with the following well known observation: 

Observation 1.1. A subgraph H of graph G has stretch factor p if and only if for every 
edge XY S G: the length of a shortest path in H from X to Y is at most p ■ \XY\. 

For three non-collinear points X, Y, Z in the plane we denote by Q)XY Z the circum- 
scribed circle of triangle AXYZ. A Delaunay triangulation of a set of points P in the plane 
is a triangulation of P in which the circumscribed circle of every triangle contains no point 
of P in its interior. It is well known that if the points in P are in general position (i.e., 
no four points in P are cocircular) then the Delaunay triangulation of P is unique [8]. In 
this paper — as in most papers in the literature — we shall assume that the points in P are 
in general position; otherwise, the input can be slightly perturbed so that this condition is 
satisfied. The Delaunay graph of P is defined as the plane graph whose point-set is P and 
whose edges are the edges of the Delaunay triangulation of P. An alternative definition 
that we end up using is: 

Definition 1.2. An edge XY is in the Delaunay graph of P if and only if there exists a 
circle through points X and Y whose interior contains no point in P. 

It is well known that the Delaunay graph of a set of points P is a spanning subgraph of 
the Euclidean graph defined on P (i.e., the complete graph on point-set P) whose stretch 
factor is bounded by C dd = 4^/3vr/9 ra 2.42 [10]. 

Given integer parameter k > 6, the Yao subgraph [16] of a plane graph G is constructed 
by performing the following Yao step at every point M of G: place k equally-separated rays 
out of M (arbitrarily defined), thus creating k closed cones of size 2n/k each, and choose 
the shortest edge in G out of M (if any) in each cone. The Yao subraph consists of edges 
in G chosen by either endpoint. Note that the degree of a point in the Yao subgraph of G 
may be unbounded. 

Two edges MX, MY incident on a point M in a graph G are said to be consecutive if 
one of the angular sectors determined by MX and MY contains no neighbors of M. 
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2. Bounded Degree Spanners of Delaunay Graphs 

Let P be a set of points in the plane and let E be the complete, Euclidean graph defined 
on point-set P. Let G be the Delaunay graph of P. This section is devoted to proving the 
following theorem: 

Theorem 2.1. For every integer k > 14, there exists a subgraph G' of G such that G' has 
maximum degree k and stretch factor 1 + 2-7r(A;cos ■ 

A linear time algorithm that computes G' from G is the key component of our proof. 
This very simple algorithm essentially performs a modified Yao step (see Section I2.3P and 
selects up to k edges out of every point of G. G' is simply the spanning subgraph of G 
consisting of edges chosen by both endpoints. 

In order to describe the modified Yao step, we must first develop a better understanding 
of the structure of the Delaunay graph G. Let CA and CB be edges incident on point C in 
G such that ZBCA <2ir/k and CA is the shortest edge within the angular sector ZBCA. 
We will show how the above theorem easily follows if, for every such pair of edges CA and 
CB: 

1. we show that there exists a path p from A to B in G of length \p\, such that: 
\CA\ + \p\ < (1 + 2vr(£;cosf)- 1 )|C5|, and 

2. we modify the standard Yao step to include the edges of this path in G' , in addition 
to including the edges picked by the standard Yao step but without increasing the 
number of edges chosen at each point beyond k. 

This will ensure that: for any edge CB £ G that is not included in G' by the modified Yao 
step, there is a path from C to B in G' , whose edges are all included in G' by the modified 
Yao step, and whose cost is at most (1 + 2ir(k cos j)~ 1 )\CB\. In the lemma below, we prove 
the existence of this path and show some properties satisfied by edges of this path; we will 
then modify the standard Yao step to include edges satisfying these properties. 

Lemma 2.2. Let k > 14 be an integer, and let CA and CB be edges in G such that 
ZBCA < 2ir/k and CA is the shortest edge in the angular sector ZBCA. There exists a 
path p : A = Mq, M%, M r = B in G such that: 

(i) \CA\ + EI^o 1 \MiM i+1 \ < (l + 27T(kco S l)' 1 )\CB\. 

(ii) There is no edge in G between any pair Mi and Mj lying in the closed region delimited 
by CA, CB and the edges of p, for any i and j satisfying < i < j — 1 < r. 

(iii) ZM^MiM^ > (^)vr, for i = 1, • • • , r - 1. 

(iv) ZCAM 1 > | -f. 

We break down the proof of the above lemma into two cases: when AABC contains no 
point of G in its interior, and when there are points of G inside AABC. We define some 
additional notation and terminology first. We define the circle (O) = QABC with center 

O, and set 6 = ZBCA. Note that ZAOB = 2G < 4vr/A;. We will use AB to denote the arc 
of (O) determined by points A and B and facing ZAOB. We will make use of the following 
easily verified Delaunay graph property: 

Proposition 2.3. If CA and CB are edges of G then the region inside (O) subtended by 
chord CA and away from B and the region inside (O) subtended by chord CB and away 
from A contain no points. 
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2.1. The Outward Path 

We consider first the case when no points of G are inside AABC. Since both CA and 
CB are edges in G and by Proposition 12.31 the region of (O) subtended by chord AB closer 
to C has no points of G in its interior. Keil and Gutwin [TO] showed that, in this case, there 
exists a path between A and B in G inside the region of (O) subtended by chord AB away 

from C, whose length is bounded by the length of AB (see Lemma 1 in |10j). We find it 
convenient to use a recursive definition of their path (for more details, we refer the reader 
to HE]): 

1. Base case: If AB G G, the path consists of edge AB. 

2. Recursive step: Otherwise, a point must reside in the region of (O) subtended 
by chord AB and away from C. Let T be such a point with the property that 
the region of Q)ATB subtended by chord AB closer to T is empty. We call T an 
intermediate point with respect to the pair of points (A,B). Let (Oi) be the circle 
passing through A and T whose center 0\ lies on segment AO and let (O2) be the 
circle passing through B and T whose center O2 lies on segment BO. Then both (Oi) 
and (0 2 ) lie inside (O), and AAO x T and ZT0 2 B are both less than ZAOB < 4ir/k. 
Moreover, the region of (Oi) subtended by chord AT that contains 0\ is empty, and 
the region of (O2) subtended by chord BT and containing O2 is empty. Therefore, 
we can recursively construct a path from A to T and a path from T to B, and then 
concatenate them to obtain a path from A to B. 

Definition 2.4. We call the path constructed above the outward path between A and B. 

Keil and Gutwin [10] , from this point on, use a purely geometric argument (with no 
use of Delaunay graph properties) to show that the length of the obtained path A = 
Mo, Mi, • • • , M r = B (where each point M p , for p = 1, ■ ■ ■ ,r — 1, is an intermediate point 

with respect to a pair (Mj, Mi), where 0<i<p<j<r)is smaller than the length of AB. 
Figure Q] illustrates an outward path between A and B. 



A = M 




Figure 1: Illustration of an outward path. 



Proposition 2.5. In every recursive step of the outward path construction described above, 
if M p is an intermediate point with respect to a pair of points (Mi,Mj), then: 

(a) there is a circle passing through C and M p that contains no point of G, and 

(b) circles Q)CMiM p and QCMjM p contain no points ofG except, possibly, in the region 
subtended by chords MiM p and M p Mj, respectively, away from C. 
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Proof. We assume, by induction, that there are circles (OmJ and (Om,) passing through 
C and Mi, and C and Mj, respectively, containing no points of G, and that the circle 
(O) = QCMiMj contains no point of G in the interior of the region R' subtended by 
chord MiMj closer to G. (This is certainly true in the base case because CA, CB E G, by 
Proposition 12.31 and by our initial assumptions). 

Since MiMj is not an edge in G, the point M p chosen in the construction is the point 
with the property that the region R of Q>MiM p Mj subtended by chord MiMj away from 
C, contains no point of G. Then the circle passing through C and M p and tangent to 
Q)MiM p Mj at M p is completely inside (OmJ U (Om.,) U R U R', and therefore devoid of 
points of G. This proves part (a). 

The region of Q)CMiM p subtended by chord MjM p and containing C is inside (OmJ U 
RU R' , and therefore contains no point of G in its interior. The same is true for the region 
of QCMjMp subtended by chord MjM p and containing C, and part (b) holds as well. ■ 

We are now ready to prove Lemma [2. 2 1 in the case when no point of G lies inside AABC. 
In this case we define the path in Lemma 12.21 to be the outward path between A and B. 

Proof of Lemma \2.S\ for the case of outward path. 

(i) With 9 = ABC A, we have | AB \ = 26 • \OA\ and sin 9 = \AB\/(2\OA\). We 

note that \CA\ + | AB | is largest when \CA\ = \CB\, i.e. when CA and CB are 
symmetrical with respect to the diameter of QCAB passing through C; this follows 
from the fact that the perimeter of a convex body is not smaller than the perimeter 
of a convex body containing it (see page 42 in [1J). If \CA\ = \CB\, sin ® = ■ 
Using elementary trigonometry, it follows from the above facts and from \CA\ < \CB\ 
that: 

\CA\ + \AB\ < \CB\ + 2e-\OA\ = \CB\ + (-^—)-\AB\ = \CB\ + (-^ 7T )-\CB\ 

sinQ cos ^ 

< (l + 2^(kcos^y l )\CB\. 
k 

The last inequality follows from G < 2ir/k and k > 2. 

(ii) If MiMj was an edge in G then, for every p between i and j, the circle Q)MiM p Mj 

would not contain C. This, however, contradicts part (a) of Proposition [231 
(Hi) If the outward path contains a single intermediate point Mi, then since M\ lies 
inside (O) = OCAB, AAM X B > vr - ZAOB/2 > vr - 2ir/k = (k - 2)n/k (note that 
ZAOB = 2 ■ ZACB), as desired. Now the statement follows by induction on the 
number of steps taken to construct the outward path between A and B, using the 
fact (proved in ^Q\) that each angle ZMj_iOjMj + i at the center of the circle (Oj) 
defining the intermediate point Mj, is bounded by ZAOB. 

(iv) This follows from the fact that ZCAM\ > ZCAB > ir/2 — ir/k. The last inequality 
is true because \CA\ < \CB\ and ZBCA < 2ir/k in ACAB. 



2.2. The Inward Path 

We consider now the case when the interior of AABC contains points of G. Let S be 
the set of points consisting of points A and B plus all the points interior to AABC (note 
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that C £ S). Let CH(S) be the points on the convex hull of S. Then CH{S) consists of 
points iVo = A and N s = B, and points N±, N s -i of G interior to AABC. We have the 
following proposition: 

Proposition 2.6. For every i = 1, • • • , s — 1 : 

(a) CNi e G, 

(b) \CNi\ < |C7JV i+1 |, and 

(c) ANi-\NiNi+i > 7r ; where ZNi-iNN+i is the angle facing point C . 

Proof. These follow follow from the following facts: CA and CB are edges in G, CA is the 
shortest edge in its cone, and hence \CA\ < \CNi\, for i = 0, • • ■ , s, and points No, ■ ■ ■ , N s 
are on CH(S) in the listed order. ■ 

Since |CiVj| < |CiVj+i| and no point of G lies inside ANCN + i (N and N + i are on 
CH(S)), CN is the shortest edge in the angular sector ZNiCN + i. Since ZNC N + i < 
ABC A < 2ir/k, by Lemma [2 . 2 1 there exists an outward path Pj between Ni and iVj+i, for ev- 
ery z = 0, 1, • • • , s— 1, satisfying all the properties of Lemma l2.21 Let A = Mo, Mi, • • • , M r = 
B be the concatenation of the paths Pj, for i = 0, • • • ,r — 1. 

Definition 2.7. We call the path ^4 = Mo, Mi, • • • , M r = B constructed above the inward 
path between A and B. 

Figure [2] illustrates an inward path between A and B. 




Figure 2: Illustration of an inward path. 

We now prove Lemma 12.21 in the case when there are points of G interior to AABC. 
In this case we define the path in Lemma 12.21 to be the inward path between A and B. 

Proof of Lemma \2.S\ for the case of inward path. 

(i) Define A" to be a point on the half-line [CA such that | C A" | = \CB\, and let 
(O") = OCA"B. Denote by a" the length of the arc of QCA'B subtended by 
chord A"B and facing ZA"CB. For every i = 0, 1, • • • , s — 1, we define arc an to 
be the arc of QCNNi+i subtended by chord iVjiVj+i and facing ZNiCNi + \. For 
every i = 0,1, ...,s — 1, we define N- to be the point on the half-line [CNi such 
that \CN[\ = |CiV i+1 |, (Oi) to be the circle OCN!N i+1 , and a- to be the arc of (d) 
subtended by chord N^N + i and facing ZN^CNi + i. Finally, for every i = 0, • • • , s — 1, 
we define N" to be the point of intersection of the half-line [CNi and circle (O"), 
and a'( to be the arc of (O") subtended by chord A r -'A r -+ 1 and facing ZN!'CN!' +1 . 
As shown in section 12.11 the length of the outward path Pj between Ni and iVj+i is 
bounded by the length of aj. Since the convex body C\ delimited by CNi, CNi + \ 
and a,i is contained inside the convex body C2 delimited by CN-, CN + i and c^, 
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by pQ, the perimeter of C\ is not larger than that of O2. Denoting by \Pi\ the length 
of path Pi, we get: 

\Pi\ < \NiN{\ + al, i = l,...,a-l. (2.1) 

Since (Oj) and (O") are concentric circles (of center C), and the radius of (Oj) is 
not larger than that of (O"), we have a- < a", for i = 0, • • • , s — 1. It follows from 
Inequality (|2.ip that: 

|*4| < lJViJV-1+oJ', i = l,---,s-l. (2.2) 

Using Inequalities (|2.ip and (|2.2h we get: 

s— 1 s— 1 s— 1 

|CA| + £ \Pi\ < \CA\ + \NiNl\ + al (2.3) 

i=0 i=0 i=0 

Noting that YX=l \ N i N i\ = \ CB \ ~ \ CA \, that Ei=o a i' = a "> and usin S the same 
argument as in part (£) of Lemma |2.2|) completes the proof. 

(ii) Since CiV p G G for p = 1, • ■ ■ , s — 1 by part (a) of Proposition 12.61 by planarity of G, 
if such an edge between two points Mj and Mj exists, then Mj and Mj must belong 
to an outward path between two points N p and N p +\ of CH{S). But this contradicts 
part (ii) of Lemma 12.21 for the case of the outward path applied to N p and N p+ i. 
(Hi) For each i = 0, • • • ,r, either Mj = Nj £ CH(S), or Mj is an intermediate point 
on the outward path between two points N p and N q in CH(S). In the former case 
ZMj_iMjM i+ i > ZNj-tMiNj+x > vr > (k - 2)ir/k for k > 14 (JVj_i and JVj- are 
points before and after M, = Nj on CH(S)), by part (c) of Proposition 12.61 In the 
latter case ZMj_iMjMj + i > (k — 2)7r/k by the proof of part (Hi) of Lemma 12.21 
applied to the outward path between N p and N q . 

(iv) This follows from \CA\ = \CM \ < \CM X \ and AACM X < ZACB < 2n/k, in triangle 
AC AM,. 



2.3. The Modified Yao Step 

We now augment the Yao step so edges forming the paths described in Lemma 12.21 are 
included in G' , in addition to the edges chosen in the standard Yao step. Lemma 12.21 savs 
that consecutive edges on such paths form moderately large angles. The modified Yao step 
will ensure that consecutive edges forming large angles are included in G' . The algorithm 
is described in Figure O 

Since the algorithm selects at most k edges incident on any point M and since only 
edges chosen by both endpoints are included in G 1 , each point has degree at most k in G' . 

Before we complete the proof of Theorem 12. 1\ we show that the running time of the 
algorithm is linear. Note first that all edges incident on point M of degree A can be mapped 
to the k cones around M in linear time in A. Then, the shortest edge in every cone can 
be found in time 0(A) (step 2. in the algorithm). Since k is a constant, selecting the £/2 
edges clockwise (or counterclockwise) from a sequence of ell < k empty cones around M 
(step 3.1.) can be done in O(A) time. Noting that the total number of edges in G is linear 
in the number of vertices completes the analysis. 

To complete the proof of Theorem 12. 1\ all we need to do is show: 
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Algorithm Modified Yao step 

Input: A Delaunay graph G; integer k > 14 
Output: A subgraph G' of G of maximum degree k 

1. define k disjoint cones of size 27r/fc around every point M in G; 

2. in every non-empty cone, select the shortest edge incident on M in this cone; 

3. for every maximal sequence of I > 1 consecutive empty cones: 

3.1. if I > 1 then select the first \£/2\ unselected 
incident edges on M clockwise from the sequence 
of empty cones and the first \£/2] unselected 
edges incident on M counterclockwise from the 
sequence of empty cones; 

3.2. else (i.e., £ = 1) let MX and MY be the incident 
edges on M clockwise and counterclockwise, 
respectively, from the empty cone; if either MX 
or MY is selected then select the other edge 

(in case it has not been selected); otherwise select 
the shorter edge between MX and MY breaking 
ties arbitrarily; 

4. G' is the spanning subgraph of G consisting of edges selected by both endpoints. 



Lemma 2.8. If edge CB is not selected by the algorithm, let CA be the shortest edge in 
the cone out of C to which CB belongs. Then the edges of the path described in Lemma \2.S\ 
are included in G' by the algorithm. 

Proof. For brevity, instead of saying that the algorithm Modified Yao Step selects an 
edge MX out of a point M, we will say that M selects edge MX. To get started, it is 
obvious that C will select edge CA. 

By part (w) of Lemma the angle ZCAM\ > vr/2 - n/k > 67r/k for k > 14. 
Therefore, at least two empty cones must fall within the sector ZCAM\ determined by the 
two consecutive edges CA and AM\, and edges AC and AM\ will both be selected by A. 
Since edge CA is also selected by point C, edge AC G G' . 

By part (iii) of Lemma 12. 2\ for every i = 1, 2, • • • , r — 1, the angle ZMj_iMjMi + i > 
{k — 2)ir/k > Wn/k for k > 12, and hence at least four cones fall within the angular sector 
ZMj_iMjMj + i. Since by part (ii) of Lemma I2.2I M;C is the only possible edge inside the 
angular sector ZMj_iMjMj_|_i, it is easy to see that regardless of the position of these four 
cones with respect to edge MjC, ends up selecting all edges MjMj_i, MjMj + i and MjC 
in steps 2 and/or 3 of the algorithm. Since we showed above that A selects edge AM±, this 
shows that all edges MjMj+i, for i = 0, • • • , r — 2, are selected by both their endpoints, and 
hence must be in G' . Moreover, edge M r _iM r = M r _\B is selected by point M r _i. 

We now argue that edge BM r _\ will be selected by B. First, observe that \BM r ^\\ < 

| AB | < |CB|. Let CD be the other consecutive edge to CB in G (other than CAf r _i). 
Because C does not select B, it follows that ZM r ^\CD < Qir/k. Otherwise, since CM r _i 
and CB are in the same cone, two empty cones would fall within the sector ZBCD 
and C would select B. Since CB is an edge in G, by the characterization of Delaunay 
edges [8], ZCM r _\B + ZCDB < tt. By considering the quadrilateral CDBM r _\, we have 



ZM r _iCD + ZDBM r ^i > tt. This, together with the fact that ZM r - X CD < Gn/k, imply 
that ZDBM r ^\ > (k — 6)ir/k > 87r/k, for k > 14. Therefore, ZDBM r -\ contains at least 



Figure 3: The modified Yao Step. 
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three cones of size 2ir/k out of B. If one of these cones falls within the angular sector 
ZCSM r _i then, since |Af r _iJB| < |CB|, BM r ~\ must have been selected out of B. 

Suppose now that ZCBM r -\ contains no cone inside and hence ZCBM r ~\ < 4ir/k. 
If one of these three cones within sector ZDBM r ~i contains edge CB, then the remaining 
two cones must fall within ZDBC and BM r —\ will get selected out of B when considering 
the sequence of at least two empty cones contained within ZCBD. Suppose now that all 
three empty cones fall within ZCBD. Then we have ZCBD > 6ir/k. 

If ZM r _\C D > 4ir/k, then since M r _iC and CB belong to the same cone, the sector 
ZBCD must contain an empty cone. Because D is exterior to Q)CBM r -.\, ZCBM r -\ < 
4ir/k, and ZM r - X CB < 2n/k, it follows that ZCDB < ZM^CB + ZCBM r ^\ < 6ir/k < 
ZDBC. Therefore, by considering the triangle ACDB, we note that \CB\ < \CD\. But 
then edge CB would have been selected by C in step 3 since the sector ZBCD contains an 
empty cone, a contradiction. 

It follows that ZM r ^CD < 4ir/k, and therefore ZM r ^BD > (k - 4)ir/k > Wir/k 
for k > 14. This means that at least four cones are contained inside sector ZDBM r -\. 
It is easy to check now that regardless of the placement of the edge BC with respect to 
these cones, edge BM r -\ is always selected out of B by the algorithm. This completes the 
proof. ■ 

Corollary 2.9. A Euclidean Minimum Spanning Tree (EMST) on P is a subgraph of C . 

Proof. It is well known that a Delaunay graph (G) contains a EMST. If an edge CB is not 
in C , then, by Lemma 12.81 a path from C to B is included in G' . All edges on this path 
are no longer than CB, so there is a EMST not including CB. m 

Since a Delaunay graph of a Euclidean graph of n points can be computed in time 
O(nlgn) [8] and has stretch factor C^ei ~ 2.42, we have the following theorem. 

Theorem 2.10. There exists an algorithm that, given a set P of n points in the plane, 
computes a plane geometric spanner of the Euclidean graph on P that contains a EMST, 
has maximum degree k, and has stretch factor (1 + 27r(fe cos f ■ Cdei, where k > 14 is an 
integer. Moreover, the algorithm runs in time 0(n\gn). 



3. Geometric Spanners of Unit Disk Graphs 

In this section we generalize our planar geometric spanner algorithm to unit disk graphs. 
Unit disk graphs model wireless ad-hoc and sensor networks and, for packet routing and 
other applications, a bounded-degree planar geometric spanner of the wireless network is 
often desired. Due to the limited computational power of the network devices and the 
requirement that the network be robust with respect to device joining and leaving the 
network, the construction/algorithm should ideally be strictly-localized: the computation 
performed at a point depends solely on the information available at the point and its d-hop 
neighbors, for some constant d (in our case d = 2). In particular, no global propagation of 
information should take place in the network. 

The results in the previous section do not carry over to unit disk graphs because not all 
Delaunay graph edges on a point-set P are unit disk edges. However, if U is the unit disk 
graph on points in P and UDel(U) is the subgraph of the Delaunay graph on P obtained 
by deleting edges of length greater than one unit, then UDel(U) is a connected, planar, 
spanning subgraph of U with stretch factor bounded by Cdei (see [Hill]). Therefore, if we 
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apply the results from the previous section to UDel(U) and observe that all edges on the 
path defined in Lemma 12.21 must be unit disk edges (given that edges CA and CB are), it 
is easy to see that Theorem 12.11 and Theorem 12. 101 carry over to unit disk graphs. The only 
problem, however, is that the construction of UDel{U) cannot be done in a strictly-localized 
manner. 

To solve this problem, Wang et al. [TT] [12] introduced a subgraph of U denoted 
LDelW(U). It was shown in [HI H] that LDeP\U) is a planar supergraph of UDel(U), 
and hence also has stretch factor bounded by C^ei- Moreover, the results in [61 [15] show 
how LDelW(U) can be computed with a strictly-localized distributed algorithm exchang- 
ing no more than 0(n) messages in total (n is the number of points in U), and having a 
local processing time of 0(Alg A) = O(nlgn) at a point of degree A. In a style similar to 
Definition [L2l LDel^ 2 \U) can be defined as follows: 

Definition 3.1. An edge XY of U is in LDel^ (U) if and only if there exists a circle 
through points X and Y whose interior contains no point of U that is a 2-hop neighbor of 
X or Y. 

We will use G = LDeffl (U) as the underlying subgraph of U to replace the Delaunay 
graph G used in the previous section. We note that G is planar, is a supergraph of UDel(U), 
and hence has stretch factor Cdei- To translate our results to unit disk graphs, we need to 
show that the inward and outward paths are still well defined in G. In particular, we need 
to show that Lemma 12.21 holds for G = LDel^ 2 \U). We outline the general approach and 
omit the details for lack of space. 

The following is equivalent to Proposition 12.31 

Lemma 3.2. If CA and CB are edges of G then the region of (O) = QABC subtended by 
chord CA and away from B and the region of (O) subtended by chord CB and away from 
A contain no points that are two hop neighbors of A, B and C. 

Proof. By symmetry it is enough to prove the lemma for the region of (O) subtended by 
chord CA and away from B. By Definition 13 .1\ there is a circle (Oca) passing through C 
and A whose interior is empty of any point within two hops of C or A. The region of (O) 
subtended by chord CA and away from B is inside this circle, so we only need to argue that 
it doesn't contain two hop neighbors of B either. If it did, say point X, then any neighbor 
of X and B would have to be a neighbor of C or A as well, a contradiction. ■ 

With this lemma in hand, the recursive construction of the outward path given in 
Subsection 12.11 can be applied to the graph G = LDel^ 2 \U). The following proposition for 
G = LDel^ (U) corresponds to Proposition 12.51 for Delaunay graphs and is proven in an 
equivalent manner: 

Proposition 3.3. In every recursive step of the outward path construction, if M p is an 
intermediate point with respect to a pair of points (Mi,Mj), then: 

(a) there is a circle passing through C and M p that contains no point of G that is a 
two-hop neighbor of C or M p , and 

(b) circles Q)CMiM p and Q>CMjM p contain no points of G that are two-hop neighbors 
of C, Mi and M p and C, Mj, and M p , respectively, except, possibly, in the region 
subtended by chords MiM p and M p Mj, respectively, away from C. 

With this proposition, we can show that Lemma 12.21 holds true for G = LDel^ 2 \U) for 
outward paths. It holds for inward paths as well, using the same argument as in Section [2721 
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Finally, it is obvious how the Modified Yao Step algorithm in Section 12.31 can be easily 
described as a strictly-localized algorithm. We can show, therefore, the following theorem: 

Theorem 3.4. There exists a distributed strictly-localized algorithm that, given a set P of 
n points in the plane, computes a plane geometric spanner of the unit disk graph on P that 
contains a EMST, has maximum degree k, and has stretch factor (1 + 2ir(k cos f • C^ei, 
for any integer k > 14. Moreover, the algorithm exchanges no more than 0(n) messages in 
total, and has a local processing time of A lg A at a point of degree A . 

Due to the strictly-localized nature of the algorithm, the algorithm is very robust to 
topological changes (such as wireless devices moving or joining or leaving the network), an 
essential property for the application of the algorithm in a wireless ad- hoc environment. 
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